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The MS charm quark mass is determined to be m c (m c ) = 1224 ± 17 ± 54 MeV from a global fit 
to inclusive B meson decay data, where the first error is experimental, and includes the uncertainty 
in a s , and the second is an estimate of theoretical uncertainties in the computation. We discuss the 
implications of the pole mass renormalon in the determination of m c . 
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I. INTRODUCTION 

Precise determinations of the bottom and charm quark 
masses are becoming increasingly important for the mea- 
surement of CKM parameters and in the search for new 
physics from experiments on B-meson and Kaon de- 
cays. For example, the theoretical uncertainty in the 
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n^vv rate is now dominated by the uncertainty 



in m c (m c ) 0. At present, the most precise determina- 
tions of the charm quark mass have been obtained from 
sum rule analyses of the e + e~ R-ratio and from lattice 
QCD (though all lattice results still use the quenched 
approximation) Ja. The results of recent analyses are 
consistent [j, |j, 0, [6j ; nevertheless it is important to con- 
sider other independent methods to extract the charm 
quark mass. Recently, fundamental standard model pa- 
rameters such as the bottom quark mass and the mixing 
angle V c b have been determined with high precision from 
inclusive semileptonic -B-decays 0,111 Hi El- In this pa- 
per, we extend the inclusive B decay analysis to also 
obtain the MS charm quark mass, m c (m c ). The charm 
quark mass was extracted previously in Ref. by the 
methods discussed here. We redo the analysis including 
a more careful consideration of renormalon effects in the 
quark pole masses. 

The experimental input consists of inclusive B de- 
cay spectra, such as the electron energy and hadronic 
invariant mass spectra in inclusive semileptonic B — > 
X c ev decay, the inclusive photon energy spectrum in 
B — > Xg'y decay, and the B meson masses and lifetimes. 
These quantities can be computed theoretically using 
heavy quark effective theory (HQET) , as an expansion in 
a s (mb) and Aqcd/to&- At present, all the quantities are 
known to order a 2 /3o> Aq CD /t7?|, and a s AQCB/nib, and 
fits to the experimental data were performed in Ref. 7] 
including all theoretical quantities to this order. We use 
the results of this fit to determine the charm quark mass. 

It is often convenient to do computations using the 
quark pole masses rn^c since they define the phase space 
boundaries in partonic perturbation theory. The heavy 
quark pole masses are known to have a renormalon am- 
biguity of order Aqcd Ujl 01 • The order Aqcd renor- 
malon means that the pole masses m& >c cannot be deter- 
mined accurately (i.e. with an uncertainty much smaller 
than Aqcd), and the a s perturbation expansion for the 



pole masses is poorly behaved, with a factorial growth in 
the coefficients at high orders. The bad behavior of the 
pole masses is already manifest in the first two orders of 
perturbation theory. In Ref. , fits to B decay spectra in 
terms of the pole mass were shown to be poorly behaved. 
Luckily, there are other quark mass definitions which do 
not suffer from a renormalon ambiguity, and which can 
be determined accurately. One of these, the lS-mass for 
the bottom quark, m^ s , was used as the bottom quark 
mass definition in the fits of Refs. 0,0- It was shown in 
Ref. [g that fits to inclusive B decays using the lS-mass 
are well-behaved, with good convergence of the perturba- 
tion series. Similarly, for the charm quark, a quark mass 
definitions without the Aqcd renormalon is the MS mass 
rn c (/i), frequently given for ii = m c , ra c (m c ). 

The difference between the bottom and charm quark 
masses is constrained by the meson mass difference, 
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where rh~B,D are the spin-averaged B, B* and D, D* me- 
son masses, and Ai,pi,ri,3 are HQET parameters. This 
equation was used to determine the bottom-charm pole 
mass difference nib — "tn c = 3.41 ± 0.01 GeV in Ref. 0, 
where we frequently suppress the superscript "pole" for 
the pole masses. The order Aqcd renormalon cancels in 
bottom-charm pole mass difference, so nib ~ ni c has no 
Aqcd renormalon, and can be accurately determined. 

The procedure for accurately determining the charm 
quark mass is simple in principle — one combines the 
values of nib — fn c and m£ s to determine m c (/i) . None of 
the three quantities suffers from a Aqcd renormalon, and 
so rn c (/x) can be determined reliably. However, it turns 
out there are some subtleties in the actual computation, 
and one must take care to ensure that the pole mass 
renormalon cancellation is properly implemented in m b — 
m c at intermediate stages of the calculation. This was 
not done correctly in Ref. pj. Including the renormalon 
cancellation properly, as done here, leads to significantly 
different results. 
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II. PERTURB ATIVE MASS RELATIONS 

The computations of the moments of the semileptonic 
decay spectra are carried out in the theory where the top 
quark is integrated out and where the bottom quark is 
treated as static. Thus the top and bottom quarks do not 
contribute in virtual quark loops, and all MS parameters 
are defined in the theory for nj = 4 running flavors. 

The 3-loop perturbative series for the relation between 
the charm pole and the charm MS mass at the scale /i m 
is IH 



TO p° lc = m {n 



g(a s , fi,m, y, m ) = 1 + e <5 (1) (a s , it, m, fi m ) 

+e 2 <5 (2) (a s ,^,m, + e 3 <S (3) (a s , /i, m, /x m ) (2) 

where a s = ai n, ~ 4 \^,) is the running coupling in the 

theory with four dynamical flavors, to = Trie 1 f ~ 4 \it, m ) is 
the MS charm mass in the theory with four dynamical 
flavors, and 5^ (a s , ^,to, fi m ) are: 
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5 (2) = ^ fl0.3193 + 2.77778 L + 

(2.98611 + 2.08333 L) L m - 0.541667 L 2 m ) , (4) 



S (3) = ^ (n6. 504 + 47.2748 i + 5.78704 L 2 + 
(6.30337 + 15.6505 L + 4.34028 L 2 ) L m 
+ (-6.20023 - 2.25694 L) L 2 m + 0.571759 L^, 



where 
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L = In A. 
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For most of the paper, the superscript (nj =4) for MS 
parameters is dropped for simplicity. The powers of the 
auxiliary parameter e = 1 indicate the respective orders 
in the loop expansion. Note that the charm quark mass 
and a s are renormalized at the scale fi m and fi, respec- 
tively. Choosing fi = [i m = m c (m c ) and a s (jn c ) = 0.39 
gives 

g = l + 0.42a s (m c )e + 1.05a s (m c ) 2 e 2 + 3.76a s (TO c ) 3 e 3 



= l + 0.17e + 0.16e 2 + 0.22e d 



(7) 



The first few terms of this series show no sign of con- 
vergence, and seem to indicate that one cannot deter- 
mine m c (m c ) to an accuracy better than, say, 20%. The 



poor perturbative behavior of Eq. J7J) is a reflection of 
the O(Aqcd) pole mass renormalon, which implies that 
the series for the pole mass behaves at high orders as 
~ fie^cts^fi^ipiPQ , where ip is the order in perturbation 
theory, /3o = 11 — 2/3n; is the one-loop coefficient of the 
QCD (3 function, n; is the number of light quark flavors, 
and /j, the renormalization scale for a s . As we will see 
explicitly in the next section, the bad behavior manifest 
in Eq. (J7J) does not limit the accuracy with which one 
can determine m c (/i) because in the relation between 
physical observables and short-distance masses, such as 
m c (/it) and mj s , the O(Aqcd) renormalon contributions 
always completely cancel. To achieve this cancellation in 
practical numerical analyses a few subtle guidelines for 
the treatment of truncated perturbative series have to be 
followed. 

For the series between the bottom quark pole and the 
bottom IS masses the higher order corrections each con- 
tain a factor of the inverse T Bohr radius ~ mj s a s in 
addition to the powers of the strong coupling. For the 
systematic cancellation of the O(Aqcd) renormalon con- 
tributions in each order of the perturbative series one 
has to use the upsilon expansion [Til ITEj . i.e. terms of 
order a™ +1 in the perturbative series between the bot- 
tom quark pole and the bottom IS masses are formally 
treated of order a™. For this reason, the powers of e and 
a s differ by one in Eq. ©. From now on, we will use e 
as the perturbation series expansion parameter, and the 
order to which we expand will be denoted by tp. 

The third order relation between the bottom pole 

and IS masses reads EE E3 («» 
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A ( 2 ) = i (n.2778 + 8.33333 Lis )+ (A (1) ) 2 (10) 

9 7T V / 
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A® = -^-(192.693+ 119.083 Lis + 52.0833 L 2 S 

9 7T V 
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The latter two terms describe the corrections arising from 
the non-zero charm quark mass. These contributions 
are included for consistency, since the charm quark is 
treated as massive in the charm pole-MS mass relation 
of Eq. Note that accounting for the charm mass cor- 
rections in Eq. © is mandatory to ensure the complete 
cancellation of the pole mass renormalon contributions in 
?7ifc — m c [l6|. This is because, in the perturbation series, 
only quarks that are treated as massless contribute to the 
bad higher order behavior that leads to the order Aqcd 
renormalon ambiguity. To be more specific, only massless 
quarks contribute to the number of light quarks ni that 
governs the asymptotic large order behavior; light quarks 
whose masses are not set to zero in the computations, as 
well as heavy quarks that are integrated out or treated 
as static, do not contribute to the Aqcd renormalon am- 
biguity in the perturbative series since their masses ef- 
fectively act as an infrared cutoff for perturbative con- 
tributions from small momenta. The charm quark has 
been treated as massive in Eq. @, and so must also be 
treated as massive in Eq. (jSJ . We discuss this issue again 
in Section |IV| The functions A^°;r L ° and A^ L ° are 
somewhat involved and can be found in [Iti, fl8| . Charm 
mass contributions first enter at order ip = 2 in the e 
expansion. We will denote order <p — 2, 3 results in- 
cluding the charm mass corrections shown in Eq. (|SJ| by 
<p = 2 C , 3 C . 



III. DETERMINING m c 

Inclusive B decay spectra have traditionally been com- 
puted in terms of the pole masses mb and m c . The for- 
mulae used in Ref. eliminated the pole mass depen- 
dence by using Eq. (Q) and Eq. ©. The resulting ex- 
pressions eliminate the Aqcd renormalon from the per- 
turbation series, and have better behaved perturbative 
expansions. The fit in Ref. used all expressions con- 
sistently to order a 2 /?o, and so used Eq. ijSJ, dropping 
the e 3 and non-BLM parts of the e 2 terms (which also 

(2) 

includes A^)- This paper focuses on the determination 
of m c (fi) , and how the results depend on the order of 
perturbation theory and the methods that are used for 
the analyses. We use the following two methods: 

Method A: Take the results of Ref. for m\ s , Ai, p\ , 

T\ and t 3 as an input, and then determine m c (m c ) using 



Method B: Take the results of Ref. for m\ s and 
m& — m c as an input, and then determine fn c (m c ) using 
Eqs. OBI- 

The difference between Methods A and B is whether 
mb — m c or the RHS side of Eq. £Q) is held fixed. Note 
that for each method the entire expressions to a given 
order in a s are kept, except for the finite charm mass 
corrections in Eq. (JSJ) as explained below. 



Note that in determining m c (/i) to order e ¥ 



0,1,2,3, we will use Eqs. © to order e v , even though 
the same expression was only used to order a 2 s (3o in the 
fit in Ref. ,?]■ Thus we break up the problem into two 
parts, determining m£ s , etc. from inclusive B decays, 
which has already been studied in detail in Ref. 0, and 
then using the results of the fit to get rn c (/j,). In both 
parts, the most precise available theoretical predictions 
are being used. Since in both parts, the leading O(Aqcd) 
renormalon contributions are being systematically can- 
celed, correlations of the behavior of perturbation theory 
in the two parts can be neglected. The errors in each part 
are determined separately, and both errors are included 
in our final error estimate. The errors in the first part of 
the problem (including the dependence on the order in 
e) were studied in detail in Refs. 00, and we use the 
same estimate here. In this work we concentrate on the 
errors in the second part of the problem. 

For method B, one substitutes Eqs. (12181) to write nib — 



fn c as 



m b - m c = f(a s )m b -g(a s )m c (n) 



(15) 



where we have shown explicitly the dependence of / and 
joiiq s , suppressing the other variables. Both f(a s ) and 
g(a s ) contain an O(Aqcd) renormalon contribution and 
have a badly behaved perturbation series. The O(Aqcd) 
renormalon cancels between the two series, and so is ab- 
sent in nib — ni c . To ensure that the series expansion 
for nib — iTic is well-behaved, it is mandatory to truncate 
f(a s ) and g(a s ) to the same order in e, when they are 
both written in terms of a s (n) at the same scale. Using 
different scales, e.g. /x = m£ s for / and \i = m c (m c ) for 
g and then truncating at the same order in e does not 
cancel the renormalon, and still leads to a badly behaved 
perturbation series which gives unreliable results. 
Solving Eq. 1|15|) for m c (/i) gives 



(16) 



f(ml s , a s (u),[i, m(p m ), Mm)™^ - (m b - m c ) 

ra ) i f^m ) 



Eqs. (HE0. 



We use Eq. (|16(l to determine m c (/i) in several differ- 
ent ways. In each case, we use Eq. (Ilfit with / and g 
expanded to order e v , with (p = 0, 1, 2, 2 C , 3, 3 C . It is im- 
portant to note that expanding to order e v means that 
all terms of higher order must be dropped. Since we 
are dealing with divergent series, one cannot retain some 
higher order terms, since they are large. The cancella- 
tions that convert the badly behaved perturbation series 
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for / and g into a better behaved series for the MS charm 
mass only take place when all terms up to a given order, 
and no terms of higher order, are included. This point is 
discussed in more detail at the end of this section. 

Our final results for m c are presented in Table [I] and 
are obtained using the following procedures (in the order 
given below): 

1. Use n = 4.2 GeV and fj, m = 4.2 GeV, to get 
to c (4.2 GeV). Eq. (|lfc>|> is expanded out analytically 
to order e v . The terms m c (/u m ) in / and g on the 
RHS are eliminated by iteration resulting in an an- 
alytic (but complicated) expression, [analytic] 

2. Use = 4.2 GeV and fi m = 4.2 GeV, to get 
m c (4.2GeV). Eq. (|TB|) is solved numerically. In 
particular, the terms In [m c (4.2 GeV)] in / and g 
are kept, [numeric] 

3. Use the result from JSJ) for m c (to c ) and determine 
m c (4.2 GeV) using the (<p + l)-loop QCD renormal- 
ization group equations. [analytic+RGE] 

4. Use the result from ijfJJl for m c (rn c ) and determine 
m c (4.2 GeV) using the (ip + l)-loop QCD renormal- 
ization group equations. [numeric+RGE] 

5. Use fi = 4.2 GeV and fi m = m c (rn c ), to get 
rn,c (rn c ). Eq. I|16(l is expanded out analytically to 
order e v . The terms m c (m c ) in / and g are elimi- 
nated iteratively, as in Q). [analytic] 

6. Use fi — 4.2 GeV and /i m = m c (m c ), to 
get m c (to c ). Eq. (|16[) is solved numerically for 
m c (m c ). In particular, the terms lnm c (m c ) in / 
and g are kept, [numeric] 

7. Use the result from JQ) for m c (4.2GeV) and deter- 
mine m c (m c ) using the (ip + l)-loop QCD renor- 
malization group equations. [analytic+RGE] 

8. Use the result from J5J for m c (4.2GeV) and deter- 
mine m c (jfi c ) using the (ip + l)-loop QCD renor- 
malization group equations. [numeric+RGE] 

For method A the same procedures are followed, ex- 
cept that one uses Eq. QJ as the starting point, and all 
results are (always) expanded to order Ag CD /m 3 in addi- 
tion to the perturbative expansion in e for the functions 
/ and g. We note that for the higher order Aqcd/™ 
corrections in Eq. JQ), we treat the bottom and charm 
mass parameter as pole masses and always eliminate the 
charm mass terms iteratively in the Aqcd/^ expansion. 
For method B, m b — m c is determined from Eq. set- 
ting the bottom mass terms in the Aqcd /nrib corrections 
equal to m£ s . Methods A differs from method B through 
the dependence of mj and m c in the 1/m 2 and 1/m 
terms in Eq. Q on the order ip. In method B, m b — m c 
is held fixed, whereas in method A, it can vary by more 
than 10 MeV between ip = and ip = 3 C . 

The reference value \i = 4.2 GeV is used in the above 
methods, since that is the value used in Ref. 0. The 



estimate of the theory uncertainty in Ref. includes half 
the a^Po term, and so already includes any uncertainties 
due to a different choice of reference scale. 



A. Renormalon Cancellation 

There are alternative ways of writing Eq. (|16|) . For ex- 
ample, one could also write the mass relations Eqs. (IL'ISI) 
in the inverse form 

rnd/j,) = gm c , ml s = fm b . (17) 

leading to four possible expressions, 

m b -m c = f ml s - gm c (fi) , (18a) 

(18b) 
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m b -m c = j ml s - g m c (fi) , 

m b -m c = / m^ s - - m c (/i) , 
9 

1 IS 1 / \ 

m b -m c — -~m b — -m c {ji) 

J Q 



(18c) 
(18d) 



If the functions /, /, g and g each are truncated to order 
(/j, only the first form represents a well-behaved series in 
practice. For the remaining three cases, it is necessary 
to reexpand the final expression in e and truncate to e v , 
leading to the first form. For the analytic procedures 
described above, all four forms lead then to the same 
answer, since the final results are always reexpanded and 
consistently truncated to e v . However, for the numeric 
procedures, where / and g are individually truncated to 
e v , and then m c is computed numerically by requiring 
that the relations analogous to Eq. (|16fl are satisfied, only 
the first form is correct. Solving Eq. (| 1 fcil) numerically 
with any of the replacements / — > 1// or g — > 1/g can 
lead to unreliable results due to badly behaved higher 
order contributions that remain uncanceled. This was 
missed in Ref. |2|; the other problem was that a s was 
not used at the same renormalization scale /i for both / 
and g. 



IV. RESULTS AND CONCLUSIONS 

The results of our analysis are given in Table [I] The 
input parameters TOj S , Ai, p\ and t\$ and their covari- 
ance matrix needed for our analysis have been obtained 
using the lS-mass fit procedure described in Ref. 0, but 
with the replacement a s (4.2GeV) = 0.22 0.2245 as 
the central value, to be consistent with a s (Mz) = 0.118. 
This makes tiny changes to the input parameters from 
the results published in Ref. Q, e.g. m£ s changes by 
6 MeV, which is much smaller than its ~ 30 MeV total 
error. The results for the 6-quark IS mass (rounded to 
10 MeV precision) remains unchanged and reads m^ s = 
4.68 ± 0.03 GeV in full agreement with the T sum rule 
analyses carried out in Refs. 0, 0] . 
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TABLE I: Results for the MS charm quark mass in MeV. The first error is statistical and includes the uncertainty in a s . The 
second error is the theory uncertainty in the fits to the inclusive spectra of semileptonic B decays carried out in Ref. 0. The 
columns are labeled by the order in the e expansion; the subscript c denotes that finite charm quark mass effects (which first 
come in at order e 2 ) have been included in the IS bottom quark mass formula in Eq. JHJ. The upper block uses method A, and 
the lower block uses method B. The different methods are discussed in the text. 
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1037 ±31 ±40 
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m c (4.2GeV) 


numeric+RGE 


1037 ±31 ±40 


893 ± 39 ± 38 


845 ± 43 ± 37 


849 ± 43 ± 37 


827 ±44 ±37 


841 ± 43 ± 37 
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analytic 


1279 ± 23 ± 39 


1210 ± 21 ± 36 


1183 ±21 ±34 


1187 ± 21 ± 34 


1172 ±20 ±33 


1188 ± 19 ±33 
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numeric 


1279 ± 23 ± 39 


1216 ± 21 ± 36 


1190 ±21 ±34 
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1192 ±19 ±33 


m c (m c ) 


analytic+RGE 


1521 ±21 ±37 


1315 ± 18 ±33 


1235 ± 17 ±32 


1240 ± 17 ± 32 


1203 ± 17 ±31 


1219 ± 17 ±31 


m c (m c ) 


numeric+RGE 


1521 ± 21 ±37 


1349 ± 19 ± 32 


1267 ± 17 ±31 


1270 ± 18 ± 31 


1226 ± 17 ±30 


1237 ± 17 ±30 



For the input parameters used in our analysis, we have 
distinguished between experimental uncertainties coming 
from the data on inclusive B decay spectra, and the the- 
ory uncertainties coming from the corresponding theoret- 
ical predictions at order a'j.Po by carrying out the fits with 
theory uncertainties (as described in Ref. 0) switched on 
and switched off. For our analysis we have also accounted 
for the uncertainty in a s (4.2 GeV) as an experimental un- 
certainty, including its correlation with ml , Ai, pi and 
Ti 3. The corresponding a s entries of the covariance ma- 
trix have been estimated from the dependence of the best 
fit for ml s , Ai, p\ and ti j3 on a s (4.2 GeV) and assuming 
that the error on a s (4.2GeV) is ±0.0114 which corre- 
sponds to an error of ±0.003 for the strong coupling at 
(i = Mz- This error has been conservatively chosen to be 
larger than the PDG estimate 0] of ±0.002 for the error 
on a s (Mz). We find that the uncertainty in a s (4.2 GeV) 
has virtually no impact on the total experimental un- 
certainty for m c (rn c ), while the experimental error for 
m c (4.2GeV) increases by 15 to 20 MeV. 

The various results shown in Table [I] are stable, and 
reflect a very good perturbative behavior in each case. 
For each entry in the table the first error is statistical ac- 
counting for the uncertainties and correlations of the in- 
put parameters m^ s , Ai, px, n,3 and a s (4.2 GeV) due to 
uncertainties in the experimental data used in the anal- 



ysis of Ref. 0, and the second error accounts for the 
theory uncertainties in the predictions used in Ref. 0. 
The results show some interesting features that are wor- 
thy of comment: It is conspicuous that the results for 
W c (jn c ) obtained directly in a single step appear more 
stable than those obtained by scaling m c (4.2 GeV) down 
to the charm mass. For example for method A, the 
two results for m c (m c ) based on the analytic determi- 
nation read (1279 - 66e - 18e^ + 6ej?) MeV [analytic] and 
(1521-203e-72e2-l5 e 3) MeV [analytic+RGE]. Both se- 
ries show very good convergence properties. On the other 
hand, the differences between the central values for the 
two series at order (l,e,e 2 c ,e 3 c ) are (242, 105,51,30) MeV. 
These values are much larger than one might guess by 
using a naive error estimate based on the good behavior 
of each individual series. This is a feature that appears 
common to the behavior of perturbation theory at the 
(relatively low) charm mass scale 4J. Since there is no 
a-priori reason to consider any single one of the series to 
be more reliable than the other, we use as an estimate of 
the theoretical uncertainty in our analysis the spread of 
central values for the different methods. 

It is instructive to compare the results at order 2, 3 and 
2 C , 3 C . While accounting for the charm mass corrections 
in the relation of the bottom pole and IS masses, Eq. © 
leads only to a small difference of at most 5 MeV at order 
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TABLE II: Linear formulas for the MS charm mass in MeV at order 3 C in terms of the input parameters for the different 
procedures used in our analysis. The input parameters are defined by Aa s = (a s (4.2 GeV) — 0.2245), Ami = (rnl s GeV" 1 — 
4.681), AA = (Ai GeV" 2 + 0.2426), Ap = ((pi -ry- r 3 ) GeV" 3 + 0.02981), 5 = {{m b ~ m c ) GeV" 1 - 3.401). The last column 
shows the maximal deviation in MeV of the linear approximation from our results over the parameter ranges |A« S < 0.015, 
|Am b | < 0.1 GeV, |AA| < 0.1, |Ap| < 0.2, |5| < 0.05. 





Method 




3„ 




max. dev. 


m„ (4 2 GeV) 


analytic 


880.7 — 


1761 Aa a + 159 AA + 865 Am b - 


- 59 Ap 


7 


rn (4 9 GpV\ 




895.4 — 


1 551 An- 4- 1 71 A A 4 858 Am,, - 


- 71 Ap 


7 


rn (4 2 GeVl 


analvtic+R(4E 


849.7 


2473 An- 4- 1 64 AA 4- 900 Ami. - 


- 61 Ap 


7 


rn (4 2 GeVl 


1± Lll±±r_*.L ±*_* |^ J. 1 V. I 1 j 


850.2 — 


2528 An- 4- 1 82 AA 4- 91 3 Ami. - 


- 76 Ap 


6 


"777 ( '! r > 1 
"t-c \ltlJCJ 


C111GL1 y U1L 


1202.0 - 


4 272 An- 4- 162 AA 4 886 Ami - 


- 60 Ap 


g 


rn c (m c ) 


numeric 


1202.6 - 


f 217 Aa a + 179 AA + 899 Am b - 


- 75 Ap 


7 


m c (m c ) 


analytic+RGE 


1232.5 - 


f 937 Aa 3 + 156 AA + 850 Am b - 


- 58 Ap 


9 


rh~ c ijn c ) 


numeric+RGE 


1247.5 + 


1126 Aa s + 168 AA + 842 Ami, - 


- 70 Ap 


8 


m c (4.2 GeV) 


analytic 




868.0 - 1920 Aa s + 836 Am b 


- 808 5 


4 


m c (4.2 GeV) 


numeric 




886.6 - 1635 Aa s + 823 Ami, 


- 798 5 


5 


m c (4.2 GeV) 


analytic+RGE 




836.6 - 2650 Aa s + 869 Am b 


- 840 5 


2 


m c (4.2 GeV) 


numeric+RGE 




840.5 - 2632 Aa s + 877 Am b 


- 850 5 


2 


m c (m c ) 


analytic 




1188.9 + 110 Aa s + 858 Am b 


- 829 5 


3 


m c (m c ) 


numeric 




1192.7 + 124 Aa s + 865 Am b 


- 838 5 


3 


m c (m c ) 


analytic+RGE 




1219.8 + 793 Aa s + 822 Am b 


- 795 5 


7 


m c (m c ) 


numeric+RGE 




1238.1 + 1052 Aa s + 809 Am 6 


- 784 5 
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2 and 2 C for m c (m c ), the difference is already as large 
as the experimental uncertainties of around 15 MeV at 
order 3 and 3 C . This behavior arises because treating 
the charm as massless in the function / of Eq. JHJ loads 
to a different 0(AQCD)-renormalon behavior in the bot- 
tom pole-lS mass series than in the charm pole-MS mass 
series where the charm mass effects are automatically in- 
cluded. Thus, as explained in Section [HI at order 2 and 

3 the O(Aqcd) renormalon contributions contained in / 
and g do not fully cancel. Therefore, to avoid this source 
of perturbative ambiguity it is mandatory to account for 
the charm mass effects in the bottom pole- IS mass rela- 
tion, and only the orders 2 C and 3 C should be considered 
for the analysis. 

To determine our final result for m c (m c ) we use the 
order results. The results obtained using method A 
and B are very similar. For the central value and the error 
estimate we therefore only use the numbers from method 
A. As the central value we use the average of the largest 
and smallest numbers. For the experimental uncertainty 
we take the largest experimental error of 17 MeV. For the 
theory uncertainty, we linearly add the error originating 
from the theory error of the input parameters ml s , Ai, 
Pi and 7+3 of 31 MeV, and the perturbative uncertainty 
of our analysis. For the latter we take half the range of 
the central values we have obtained, which is 23 MeV. 
This gives 

m c (m c ) = 1224 ± 17 cxp ± 31 B ± 23 mc MeV , (19) 



with the subscripts denoting the sources of the various 
uncertainties. We add the two theoretical uncertainties 
linearly, to be conservative. Moreover, this more conser- 
vative treatment should also account for the uncertainty 
arising from the unknown order a s Aq CD /m corrections 
in Eq. The latter uncertainty is reflected by the vari- 
ation of nib — n^c with the order ip, and is about 10 MeV. 
Thus our final result is 

m c (m c ) = 1224 ± 17 ± 54 MeV. (20) 

The first error is the experimental 1-cr error and the sec- 
ond is the theoretical error. Note that the theoretical 
error should not be interpreted as a 1-a error, since it is 
not statistical. Rather it represents the range in which 
we believe that the true charm mass is located with high 
probability (much greater than the 67% probability of 
a 1-er interval). Our result is consistent with the most 
recent result for m c (rn c ) from low-n e + e" moment sum 
rules Q. However, our result is somewhat lower than 
the most recent results obtained from lattice QCD 0, 0| , 
which were obtained in the quenched approximation. 

Linear approximation formulae showing the depen- 
dence of our results on the input parameters m^ s , Ai, 
pi, 7+3, a s (4.2 GeV) and mb — m c for methods A and B 
at order ej? are displayed in Table ITT1 

In our analysis we have also determined m c (4.2GeV) 
using the same procedure as used for m c (m c ), 

m c (4.2GeV) = 873 ± 41 oxp ± 34 B ± 23 mc MeV .(21) 
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Our final result is 

m c (4.2GeV) = 873 ± 41 ± 57 MeV. (22) 

While the theoretical uncertainties are comparable to the 
corresponding uncertainties for m c (m c ), we find that the 
experimental errors are 15 to 20 MeV larger. As men- 
tioned before, this increase originates from the uncertain- 
ties in a s which have a rather large impact on the un- 
certainty in m c (4.2GeV), but almost none for m c (m c ). 
The origin of this behavior can be clearly seen from the 
linear approximation formulae shown in TablelTTl the val- 
ues of m c (m c ) and m c (4.2 GeV) both increase for larger 
input values m£ s ; on the other hand, while the results 
for m c (m c ) increase with a s , the values m c (4.2GeV) 
decrease. The fits of the semileptonic inclusive B de- 
cay spectra show that to^ s becomes smaller for larger a s 
values, and so the Atoj, term tends to cancel the Aa s 
term for m c (m c ), but comes in with the same sign for 
m c (4.2GeV). 

Based on our error analysis it is straightforward to 
discuss how the uncertainties in our charm mass anal- 
ysis might be further reduced in the future. More pre- 
cise experimental data on the semileptonic inclusive B 
decay spectra can not reduce the error significantly be- 
cause they constitute only a relative small contribution 
to the total error. In this respect our analysis is similar 
to sum rule analyses based on the e + e~ R-ratio and lat- 



tice QCD where theory uncertainties dominate. A more 
precise input value for a s will not affect the uncertainties 
of m c (m c ), but can reduce substantially the experimen- 
tal error of m c (4.2GeV). A full 0(a 2 s ) analysis of the 
semileptonic inclusive B decay spectra would reduce the 
theory uncertainty of the input parameters ml s , Ai, p\ 
and Ti ,3 of our charm mass analysis. Moreover, a deter- 
mination of C(a s AQ CD /m 2 ) corrections would reduce the 
scheme dependence of the Aq CD /to contributions in the 
OPE for the B-D meson mass difference. We believe that 
these computations are achieveable with present technol- 
ogy and could reduce the theory errors shown in Table |U 
by at least a factor of two. Finally, a computation of 
the full order corrections to the bottom-charm pole 
mass difference as a function of the MS charm mass and 
the bottom IS mass could further reduce the theoretical 
uncertainty in Eq. I|l(j|) . Such work would allow one to 
reduce the uncertainty of rn c (rn c ), with a conservative 
estimation of theory errors, to well below the 50 MeV 
level. 
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